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ABSTRACT 

The Space Interferometry Mission (SIM) will observe sources in crowded fields. Recent work has 
shown that source crowding can induce significant positional errors in SIM's astrometric measurements, 
even for targets many magnitudes brighter than all other crowding sources. Here we investigate whether 
the spectral decomposition of the fringe pattern may be used to disentangle the overlapping fringes 
from multiple blended sources, effectively by performing synthesis imaging with two baselines. We 
find that spectrally dispersed fringes enable SIM to identify and localize a limited number of field 
sources quite robustly, thereby removing their effect from SIM astrometry and reducing astrometry 
errors to near photon noise levels. We simulate SIM measurements of the LMC, and show that (a) SIM 
astrometry will not be corrupted by blending and (b) extremely precise imaging of mildly crowded fields 
may be performed using only two orthogonal baseline orientations, allowing microarcsecond positional 
measurements. We lastly illustrate the method's potential with the example of astrometric microlensing, 
showing that SIM's mass and distance measurements of lenses will be untainted by crowding. 



1. INTRODUCTION 

The Space Interferometry Mission (SIM) will perform 
astrometry with unprecedented precision, measuring posi- 
tions at the microarcsecond level. This will allow numer- 

' ous astrophysical experiments to be performed that pre- 

, viously were impossible. Even though its 10 meter base- 
line allows resolution only of ~ 10 mas, SIM achieves mi- 

' croarcsecond astrometric precision by centroiding fringes 
to better than 1 part in a thousand. If SIM can meet 
its technical specifications, the primary source of error for 
single star astrometry would be photon noise, which can 
be beat down as 1/VN until the systematic error floor is 
reached. For multiple stars in the field of view, however, 
other errors may be incurred. For example, the extra flux 

1 of photons from field sources gives additional photon noise, 
which increases the amount of on-source integration time 

i to achieve a given astrometric precision. Another source of 
error, which cannot be reduced by further integration time, 
is the distortion in the measured fringe pattern caused 
by the superposition of fringes from multiple incoherent 
sources. Since SIM achieves its revolutionary precision 
by fitting the observed fringe pattern in exquisite detail, 
even small scale deviations from the single star fringe pat- 
tern can lead to large astrometric errors. Several possible 
SIM experiments, such as microlensing, globular cluster 
dynamics, or stellar motions in external galaxies like M31, 
require SIM astrometry in crowded fields, and so it is nec- 
essary to determine whether or not crowding will corrupt 
SIM measurements irretrievably. 

One might imagine that SIM astrometry of targets many 
magnitudes brighter than other sources in the field of view 
would be unaffected by the presence of those sources. How- 
ever, Rajagopal et al. (2001) have recently shown that 
source confusion in crowded fields can lead to significant 
errors (~ few /Lias) in individual position measurements, 
even when the confusing field sources are just a few per- 
cent of the target's brightness. Such errors are comparable 



to those expected from photon noise and systematic error, 
and could pose a serious threat to astrometric measure- 
ments in crowded fields if uncorrected. Brighter sources, 
such as binary companions or highly luminous field stars 
> 5% in brightness, could completely destroy SIM astrom- 
etry. To correct these measurements, one would like to 
know the positions and brightnesses of all the important 
sources in the field of view - that is, one would like an 
image. A possibility is to obtain, prior to SIM observa- 
tions, either space telescope or ground based adaptive op- 
tics (AO) snapshots of the field. Unfortunately, as we show 
below, the rough positions obtainable (~ 40 mas) from 
such images are of insufficient precision to rectify SIM 
astrometry. Another alternative is to perform synthesis 
imaging with SIM. In analogy to radio interferometry, one 
might orient SIM along numerous baselines, using multi- 
ple siderostat pairs to fill in the (u, v) plane (Allen 2000). 
Note, however, that a possible SIM design currently under 
consideration involves just one 10 meter baseline. A third 
alternative, which we discuss here, is to use SIM in its as- 
trometry mode. If the sources all lie several fringes away 
from each other then the (it, v) coverage obtained from 
just two baselines, using multiple frequency channels, is 
sufficient to disentangle the sources from each other and 
measure all of their positions to few /ias precision. This is 
a marked increase in astrometric precision relative to that 
possible from SIM imaging, ~ 250 /ias (Allen 2000), and 
thus allows significantly improved proper motion measure- 
ments. We show that, for the currently expected charac- 
teristics of the SIM satellite, and for typical levels of crowd- 
ing in the LMC, Bulge, and the outer regions of globular 
clusters, these two orthogonal baselines are adequate to 
achieve excellent astrometric or proper motion precision. 
Note that similar ideas have been used previously in sim- 
ulations of closure phase measurements (Schloerb 1990) 
and in visibility measurements of the binary star Capella 
(Koechlin et al. 1979). 
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To see why just two baselines are sufficient, consider a 
synthesis imaging measurement of several point sources. 
This involves measuring the complex visibility V at nu- 
merous points in the (u,v) plane, where u = B x /X, 
v = By/X are spatial frequencies, and B x and B y are 
the projections of the baseline separation along the x, y 
axes respectively. The visibility V(u, v) and brightness 
distribution F(x,y) are related by Fourier transforma- 
tion, e.g. V(u, v) — J dxdyF(x, y) exp(— 2iri[ux + vy]), via 
the van Cittert-Zernike theorem (see for example Lawson 
2000, Chapter 2). For an arbitrary brightness distribu- 
tion F(x,y) sufficiently many u,v points must be sam- 
pled to permit accurate Fourier inversion of the visibil- 
ity. A typical astrometry measurement relies upon just 
two baseline orientations, which would appear insufficient 
to localize more than one source position (x and y co- 
ordinates). Note, however, that each frequency channel 
at a given baseline orientation constitutes an indepen- 
dent u, v measurement, since u = kB x /2ir, v = kB y /2ir 
(Schlocrb 1990). The multifrcquency u, v coverage from 
two orthogonal SIM baselines is limited to points on the 
u, v axes (see Figure 1), and is clearly insufficient to re- 
construct an arbitrary brightness distribution. For point 
sources, on the other hand, this coverage suffices. Consider 
F{x,y) = J2i Fi(k)5(x - Xi)S(y - yi) corresponding to a 
sum of point sources, each with spectrum Fj (k) . Taking for 
now each Fi independent of k gives the measured visibil- 
ity as V(u, v) — Fi exp(— 2-ni[xiU + y-iv]). Now, instead 
of the usual 2D Fourier transform to recover the bright- 
ness distribution, consider a single ID Fourier transform, 
say along the u direction. This gives fi(v)5(x — x{), 
and similarly for the ID transform along the v direction. 
Clearly, the sources' x coordinates are measurable just by 
sampling the visibility at multiple u along a line of con- 
stant v, and similarly for the y positions. By placing each 
source on a color-magnitude diagram (CMD) using the 
Fi(k), the x sources can be matched up with the y sources, 
giving the unique 2D coordinates of each bright source in 
the field (Dyck et al. 1995). If the pairing is not obvious 
from the CMD, then a third baseline orientation is required 
to match up the x and y sources, but note that this need 
be done only once, at any time. In theory, therefore, the 
two orthogonal SIM baselines required to perform single 
source astrometry are all that are necessary to perform 
imaging of point sources. In practice, finite bandwidth, 
Nyquist sampling and photon noise will limit point source 
resolution. We explore these effects in the next section. 
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Fig. 1. — Schematic depictions of (u,v) coverage for (a) synthe- 
sis imaging with many baseline orientations, (b) astrometry with 
two baseline orientations, and (c) astrometry with two baselines and 
channeled spectrum. Only half of the (u, v) plane is illustrated since 
the brightness distribution is real; redundant points are not plotted. 



2. FITTING DISPERSED FRINGES 

The above procedure of Fourier transforming along the 
u or v direction is equivalent to taking the spectrally 
dispersed ID fringe pattern for each orthogonal base- 
line orientation, and Fourier transforming with respect 
to frequency k. This is the same procedure that is of- 
ten used for single source phase and group delay esti- 
mation (Lawson 2000, Chapter 8). Let F(k) equal the 
product of source spectrum, throughput, and detector re- 
sponse. Note that since CCD's count numbers of pho- 
tons, not power, F has units of numbers. The interfer- 
ometer response at frequency k to a source offset from 
the phase center by an angle 9 along the baseline B is 

^-^(1 + cos kd) — ^-lr^(l + cos kBO) where d is the phys- 
ical delay. Fourier transforming with respect to k gives 
[2F(l) + F(l -d) + F(l + d)]/4, with I the Fourier conju- 
gate to k and F(l) the Fourier transform of F(k). In most 
cases, F(k) is slowly varying, and so its Fourier trans- 
form has power mainly at low frequencies. For example, 
for F(fc)=const, F(l) = 5(1) and the Fourier transform of 
the fringe pattern becomes [26(1) + 5(1 — d) + 5(1 + d)]/4. 
Likewise, for a rectangular bandpass, we would have a sine 
function. For simplicity, we will focus on rectangular F(k), 
although the generalization is clear. 

As expected, each source in the field of view appears 
as a narrow peak in the Fourier transform space (actually 
a pair since the brightness is a real quantity). Also note 
that half of the power lies at zero frequency, due to the 
1/2 which offsets the average fringe intensity from zero. 
Figure 2 illustrates the pattern for two point sources. 
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Fig. 2. — Fourier transform of the channeled spectrum with two 
equally bright sources, one at the phase center and the other at 250 
mas. The abscissa is the Fourier conjugate to wavenumber, the or- 
dinate is delay, and the numbers along the axes label pixels. Note 
that sources appear as lines in the Fourier transform, which cross 
at the position (delay) of the respective source. Additional sources 
should appear as additional, distinct lines. This figure neglects pho- 
ton noise, which can hide the signal from faint sources. 
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Now, assume that SIM uses n c frequency channels. For 
simplicity, we take each channel to have the same fre- 
quency width. This assumption allows the use of the 
FFT, otherwise a more general discrete Fourier transform 
or Lomb-Scargle periodogram would be used. With this 
assumption, the Nyquist theorem tells us that each bin in 
Fourier space has a size equaling the white-light coherence 
envelope Siz6 A max A m i n / A m ax A m i n )1 which for a rect- 
angular 0.4 — 1.0 fim bandpass and 10 meter baseline is 
roughly 13.75 mas. This makes sense, since SIM should 
be able to resolve sources separated by more than one 
coherence envelope size. Clearly, decreasing the effective 
bandwidth should diminish SIM's ability to resolve point 
sources. Varying the number of channels should affect 
only aliasing. With fewer, broader channels, a source near 
the edge of the field of view has reduced fringe contrast 
and would be mistaken for the combination of a fainter 
source near the middle of the FOV and uniform surface 
brightness. This could corrupt astrometric measurements 
of sources towards the edge of the field of view. 

We have simulated SIM astrometric measurements of 
crowded fields with the above considerations. We assume 
that SIM scans repeatedly over a small range in delay, cor- 
responding to the central fringe size at the mean frequency, 
which for a 10 m baseline and 0.4 — 1.0 /im bandpass is 
11.8 mas. We divide this range into rid bins 1 mas large, 
and integrate the fringe pattern f(k,d) = F, (1 + cos kd)/2 
over rid x n c total bins, where n c = 64 is the number of 
wavelength channels spanning the bandpass. We sum the 
contributions for all sources in the field of view, and ap- 
proximate photon noise by adding to each bin a random 
number of photons drawn from a Gaussian with a = >/~N. 
We neglect all instrumental noise, such as read noise, dark 
current, thermal drift, etc., basically assuming that SIM 
is a perfect astrometer. These additional errors may be- 
come important for the faintest sources. To simulate SIM 
astrometry and photometry, we fit the resulting fringe pat- 
tern to the form 

l + coskB(6-6i) 



)dd 



= N - 



N 4 - 



The number of sources, their positions and brightnesses 
are first estimated using the Fourier transform method 
described above. Starting with these initial values we fit 
by minimizing the error function y 2 = J2i,j( N ij, data ~ 
model) 2 /Nij. data- The fitted parameters are each 
source's brightness Ni and position 8i, as well as uniform 
surface brightness N . This assumes that the product of 
the source spectrum, throughput, and detector response is 
constant with respect to wavelength, but it is not difficult 
to generalize to include simple spectral dependence. 

An estimate of the centroiding precision possible by such 
fringe fitting may be computed using the minimum vari- 
ance bound, hereafter MVB (Gould 1995). When fitting 
data to a function f(x) with normalization and one ad- 
ditional parameter x , the minimum variance expected in 
the fitted xo is 

-l 



var(xo) 



1 

N 




din/ 

dxn 



where (g) = (J fgdx)/(J fdx). Using the MVB formula, 
and taking TV total photons from a single point source, we 



expect centroiding errors ere = 1-87 iV -1 / 2 mas for scans 
covering 11 mas and the wavelength range 0.4 — 1 /xm. 
Using the white- light fringe pattern (n c —l) we expect 
o B = 2.91V- 1 / 2 mas. So for V = 340,000 photons we ex- 
pect 3.2 /zas and 5 /j,as errors for the spectrally dispersed 
and white-light fringe pattern, respectively. This analytic 
estimate compares well with a Monte Carlo calculation of 
the expected astrometric errors using the above fringe fit- 
ting procedure. Using bins of size 1 mas (rid ~ H) gives 
errors of 3.8 /Ltas and 5 /Ltas, in agreement with the MVB 
estimate for white-light fringes and somewhat larger than 
the MVB estimate for the spectrally dispersed fringe pat- 
tern. Increasing the number of bins to ~ 100 gives Monte 
Carlo estimates in agreement with the MVB estimate. One 
might conclude from this that it is always better to cen- 
troid using spectrally dispersed fringes rather than broad- 
band, white-light fringes since the errors are smaller for a 
given number of photons. However, keep in mind that we 
have neglected additional errors such as read noise, and 
there may be a reduction in throughput when the fringes 
are spectrally dispersed. 

3. FITTING MULTIPLE SOURCES 

Having confirmed that our Monte Carlo errors agree 
with analytic estimates for single sources, we now explore 
errors for multiple sources. As mentioned earlier, addi- 
tional sources in the field of view can increase astrometry 
errors by contributing extra photon noise, and by distort- 
ing the fringe shape. Fitting for multiple sources cannot 
remove photon noise, but it can mitigate the distortion 
arising from superimposed fringes. To ease the compari- 
son with unblended centroiding as described in the pre- 
vious section, we set the on-source integration time to 
equal the time necessary to accumulate: 340,000 photons 
from the science target. In terms of a signal-to- noise ratio 
(S/N), this holds fixed the signal while letting the noise 
vary. When gauging how well fitting can remove the ef- 
fects of blending, we should compare to a case with the 
same ratio of signal (photons from science target) to noise 
(square root of total number of photons). For example, 
consider centroiding a single star in the presence of uni- 
form surface brightness that contributes an equal number 
of photons as the star. Since the S/N is smaller by a fac- 
tor of y/2 relative to the isolated star case, the centroiding 
errors will be greater by a factor of 

As our first example, we consider measurements of a 
bright microlensing event in the LMC. We place a V = 19 
target at the phase center, and randomly place field stars 
in the 800 mas FOV. The number of field stars, and their 
luminosity function (LF), are chosen to match the crowd- 
ing and LF of the LMC derived by the MACHO project 
(Alcock et al. 2001). On average, MACHO finds 10-16 
field stars within 1-2" of each observed star. To be conser- 
vative, we have used 16 field stars in all of our simulations. 
The average brightness of a star drawn from this LF is ap- 
proximately 1.4% the brightness of a V — 19 star, so by 
adding 16 of them, we add approximately 22% more pho- 
tons. We thus expect a roughly 10% increase in astrom- 
etry errors from photon noise alone. Since we normalize 
the on-source observation time to give us 3.8 /xas photon 
noise errors from the V = 19 source alone, we expect 4.2 
/j,&s errors from photon noise. Any additional errors will 
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Fig. 3. — Histogram showing position residuals for multiple source fitting to V = 19 source and 16 stars drawn from LMC luminosity 
function. The first panel shows the histogram for random placement of stars. The peaks at the edge of the plot are due to events falling 
outside the histogram range. The second panel shows the histogram for the same calculation, restricting sources to lie more than 25 mas from 
the target at the phase center. See the text for more details. 



be due to source confusion (Rajagopal et al. 2001). 

We follow the same fitting procedure described above, 
allowing for multiple sources at different positions and 
with different brightnesses. Note, however, that faint 
sources may not be detectable above photon noise, so we 
will only be able to remove the brightest sources. Since 
in this example we are only interested in the measure- 
ment errors for our single target, we only report measure- 
ment errors for the V = 19 target. Panel (a) of Figure 
3 shows the histogram of position residuals for this sim- 
ulation. Overall, the standard deviation of this error dis- 
tribution is ag = 10.5 /ias, much larger than the 4.2 /ias 
errors expected from photon noise. Note, however, that 
the histogram may not be well described by a single "av- 
erage" error value. We find that the error distribution is 
characterized by a narrow core superimposed on a slowly 
decaying tail of outliers. This histogram is quite similar to 
that found by Rajagopal et al. (2001), although for differ- 
ent reasons. In their simulations of white- light astrometry 
in crowded fields in M31, the LMC and Galactic Bulge, 
the long tail was due to the occasional bright field star. 
Here, in contrast, bright field stars are easily identified 
and removed using dispersed fringes. The long tail here is 
due to the occasional star landing nearby our target and 
becoming unresolved from it. To demonstrate this point, 
panel (b) of Figure 3 shows the histogram for the same 
simulation, but here we restrict field stars to lie more than 
25 mas from our target. Note that the long tails have van- 
ished, and that the distribution's standard deviation has 
significantly diminished, to 5.3 /zas, near to the value ex- 
pected from photon noise. Clearly, almost all of the 10.5 
/ias errors in this case arise from unresolved or marginally 
unresolved stars falling close to the target. Rajagopal et 
al. (2001) have shown that these quite significant errors 
largely cancel for proper motion measurements, and so for 
many SIM applications they are not a worry. For mi- 
crolensing, however, confusion errors will be a problem in 
crowded fields. This is because microlensed stars undergo 
different proper motions than other stars in the field (even 



binary companions), and their brightnesses vary during 
the course of the lensing event. Therefore, the cancellation 
of errors that occurs for typical proper motion measure- 
ments will not occur here, and we expect proper motion 
errors of the same magnitude as the large position errors 
found above. Fortunately, as we show below, it is still pos- 
sible to disentangle the blend effect using the complete set 
of position measurements. Thus, even in this most badly 
affected case, it is still possible to recover accurate masses 
and distances from astrometric microlensing events. 

Generally, blend stars in SIM observations will be far 
fainter than the intended science targets, which are chosen 
for their brightness. However, it is still worthwhile to con- 
sider what effects comparably bright sources would have. 
To quantify this, we simulate SIM measurements of two 
equally bright stars in the field of view. Again we normal- 
ize the on-source integration time to give 340,000 photons 
from the science target. Photon noise alone should thus 
give us errors of -s/2-3.8 /ias = 5.4 /ias. We plot the results 
of our simulations in Figure 4. We can distinguish three 
kinds of blends from this figure. First, note that bright 
blends falling outside the central fringe are resolved by 
SIM and can be removed. The offsets due to such sources 
are consistent with zero, and the dispersions (~ 6 /ias) are 
not much larger than what is expected from photon noise 
alone. Clearly, such blends pose no threat to SIM astrom- 
etry using spectrally dispersed fringes. The average offset 
is not zero, however, for sources falling within the central 
fringe. Sources at distances < 2 mas are completely unre- 
solved by SIM. As expected, then, the measured centroid 
falls at the center of light, which is 0.5 mas for an equally 
bright blend at 1 mas, and 1 mas for an equally bright 
blend at 2 mas. Note that although the offset is large, 
the standard deviation is small. This is just what we ex- 
pect; two equally bright sources at and 1 mas should 
appear to SIM to be a source twice as bright located at 
0.5 mas, and so the errors should be smaller by a factor 
of 2 -1 / 2 compared to single source astrometry. The only 
blend positions giving significant errors are those that are 
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marginally unresolved, i.e. falling barely inside the central 
fringe, at separations of 3-4 mas. There would seem to be 
no way to handle such blends with a single baseline ori- 
entation, however an additional pointing along a different 
baseline orientation could resolve the two sources. 

ft is apparent from Figure 4 that the vast majority of 
bright field stars will pose no threat to SIM astrometry. 
The fact that the fit residuals are so small, comparable to 
photon noise, means that fitting can reliably remove the 
fringes from blending sources in the field of view. Since 
we are able to recover /xas astrometry even in the pres- 
ence of blends, this must mean that we are able to fit the 
blend positions to microarcseconds themselves. Since we 
are doing this for all important blends in the field of view, 
we are basically imaging the field of view. It is therefore 
interesting to explore the imaging capabilities of SIM us- 
ing just the spectrally dispersed fringes measured along 
two orthogonal baselines. From the discussion above, we 
expect that we can get reliable positions as long as the 
sources are separated by more than one white-light coher- 
ence envelope. In terms of the Fourier transform space, we 
expect to measure reliable positions for sources separated 
by more than one pixel. With n c wavelength channels, we 
have n c /2 distinct pixels in the Fourier conjugate variable 
since the brightness is real. We have assumed 64 chan- 
nels, so we have 32 pixels. Assuming that sources should 
be separated by 3 pixels for reliable positions means that 
we expect that for up to 10 sources in the FOV, we can 
measure positions reliably. To test this, we simulated SIM 
astrometric measurements of 10 stars in the field of view. 
We first simulated 10 stars of equal brightness; the results 
are shown in panel (a) of Figure 5. We find that the as- 
trometric errors scale like t~ 1 / 2 where t is the on-source 
integration time, so they are mainly due to photon noise. 
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Fig. 4. — Fit residuals for astrometric measurements of a single 
star, with an equally bright blend in the field of view, as a function 
of distance to the blend. The thin curve is the average offset of the 
fit source position from the input position and the thick curve is the 
standard deviation of those offsets. 



For the simulation plotted in panel (a), the astrometric 
errors were a — 6 ^as. In panel (b), we plot the scaling 
of errors with brightness. Using the above arguments on 
S/N, we expect the errors to scale inversely with each star's 
brightness, and this is verified in Figure 5(b). Clearly, far 
smaller errors are attainable using dispersed fringes than 
the ~ 250 /uas errors expected from synthesis imaging of 
many baseline orientations (although keep in mind that 
we have neglected errors such as thermal drift). Thus, it 
is possible that highly precise astrometric imaging may be 
achieved using only two SIM pointings. 

4. APPLICATION: astrometric microlensing 

One possible application of SIM astrometry in crowded 
fields is microlensing. There arc numerous astrophysical 
applications of SIM observations of microlensing events 
(Paczyhski 1998; Boden et al. 1998; Safizadeh et al. 1999; 
Gould & Salim 1999; Salim & Gould 2000), most of which 
entail crowded fields. For ordinary proper motion, Ra- 
jagopal et al. (2001) have shown that astrometric errors 
arising from source confusion largely cancel. This is be- 
cause any offset caused by additional sources remains ba- 
sically constant as the target executes its proper motion, 
and cancels when the relative motion of the target is mea- 
sured. However, this is not true for microlensing. The 
reason for this is that microlensing changes the target's ap- 
parent brightness while also causing the apparent proper 
motion. The astrometric offset caused by blending, there- 
fore, does not remain constant and will not cancel. Thus, 
we might expect that blending could seriously damage mi- 
crolensing measurements, rendering SIM observations use- 
less in crowded fields. Fortunately, the results discussed 
in previous sections indicate that blending by bright stars 
should not corrupt SIM astrometry, except when those 
sources fall within the white-light coherence envelope of 
the target star. However, binary companions to lensed 
stars could easily fall inside the central fringe. Since fringe 
fitting could never detect such stars, just as fringe fitting 
cannot resolve the two microlensed images, our concerns 
about crowding could still hold. Fortunately, as we de- 
scribe below, even sources falling within the central fringe 
may be detected and removed, not from individual SIM 
measurements, but by using the entire set of observations 
of the microlensing event. 

As discussed in the previous section, field sources may 
be detected and essentially removed as long as they lie 
outside the central fringe of the target. In fact, even 
marginally resolved sources falling barely within the cen- 
tral fringe may be detected and removed if bright enough, 
by the distortion they cause in the fringe shape. The only 
sources undetectable are those which do not distort the 
fringe shape, but merely shift its centroid. These sources, 
lying well inside the central fringe, shift the fringe cen- 
troid to the center of light of the target and blend, and 
add all of their light to the target. What SIM would de- 
tect would be a star of brightness equaling the sum of the 
target and blend's brightness, located at their center of 
light. For this reason, multiple blends falling inside the 
central fringe arc equivalent to a single blend, and so we 
consider only the single blend case. Such blends affect the 
microlensing lightcurve and astrometric motion in a simple 
way (Han & Kim 1999), and we show that SIM observa- 
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astrometric error [mas] relative source brightness 

Fig. 5. — Simulated astrometric imaging errors for 10 sources. Panel (a) show the error histogram from a simulation of imaging of 10 
equally bright sources in the 800 mas field of view. The standard deviation of this histogram is a = 6 ^tas. Panel (b) shows the scaling of 
error with brightness for a simulation of 10 sources with varying brightnesses. Plotted points are errors as a function of brightness relative to 
the brightest star in the field. The dashed line is not a fit, but shows the expected scaling: error inversely proportional to brightness. 



tions may be corrected without degeneracy to extract the 
desired measurables. 

It is first necessary to review astrometric microlensing; 
see Boden et al. (1998) for more details. Microlensing of 
a single source by a single lens gives two images of the 
source. If the source is located at position x s = 9^u rela- 
tive to the lens, then the two images appear at positions 
with magnifications 



£1,2 



i{u± \/m 2 + 4)/2 



2uyv 



1 

1** 



Here, 9^ is the angular Einstein radius, u = \u\ is the lens- 
source separation in units of 9e (not the u from the u, v 
plane), and positions are relative to the lens position. SIM 
cannot resolve these images, so in the absence of blending 
it detects a source with position x — x s + Ax and magni- 
fication A given by 



Ax = 



A = 



ly/u 2 + 4 



Note that Ax traces out an ellipse. 

It is now easy to see how unresolved blends affect this. 
If the lensed source has unlensed brightness F s , and the 
blend has brightness F b and position x b , then SIM detects 
a source at position i* Q bs with brightness F to t given by 



AF S + F b , 



*^obs 



AF s (x s + Ax) + F b x b 



An example of this is given in Figure 6. From the illustra- 
tion, it may appear that blending would ruin SIM astrom- 
etry, since the observed motion bears little resemblance to 
the unblended motion. However, our simulations indicate 
that it will be possible to recover the unblended motion. 
This is largely because SIM's exquisite photometry (recall 
N > 340, 000) permits the nondegenerate decomposition 



of the photometric lightcurve into lensed and unlensed 
(blend) light. That is, F s , F b , and A may be measured 
from the photometry alone. Given the blend fraction, we 
now show that it is easy to fit the astrometric data to 
recover the microlensing parameters 9e, £e = | | 1 , etc., 
which are then used to determine the lens mass and dis- 
tance. 
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Fig. 6. — Microlensing source and motion of lens and image 
center-of-light. The large dot shows the location of the blend star, 
the small dot the location of the target star. The straight arrow 
shows the location and relative motion of the lens near the peak of 
the event, with x's placed every week, = 20 days and 6>e = 0.32 
mas. The elliptical curve asymptoting to the target source shows the 
unblended astrometric motion of the image center of light, while the 
irregular curve is the observed motion of the center of light including 
the blend star light. Arrows show the direction of motion. 
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Fig. 7. — Fit residuals for microlensing event. Each point in the 
scatter plot represents a different simulated event; plotted are the 
residuals in fitted 6*e and tE- 10,000 simulations of unblended mi- 
crolensing events are plotted in the first panel; the standard de- 
viations are ag = 0.009 mas and at = 0.1 days, for input values 
0E = 0.32 mas and <e = 20 days. In the second panel, we plot 1000 
simulations of microlensing events with a V = 19 target near the 
phase center and 16 additional stars drawn from LMC luminosity 
function placed randomly in FOV. Here, the standard deviations are 
erg = 0.008 mas and at = 0.14 days for the same input parameters 
as the unblended case. The last panel plots errors for 10,000 simula- 
tions of microlensing events with 16 field stars, but with white-light 
fringe fitting. Here, $e is systematically underestimated by 0.03 mas 
(roughly 10%), with a dispersion of 0.017 mas, while tE has a dis- 
persion of 0.44 days. 



The blended microlensing event is parametrized by 11 
quantities: the above mentioned #e, £e, x s , Xb, F s , and Fb, 
as well as the normalized impact parameter u m i n , time of 
closest approach to: and relative proper motion direction 
<p. For simplicity we assume that the same blend bright- 
ness Ft applies for both the x and y astrometric measure- 
ments. This will not be true in general (a field star whose 
projected position relative to the target happens to fall 
inside the central fringe under one orientation need not 
fall inside the central fringe along the orthogonal orienta- 
tion). However, it is trivial to generalize our treatment to 
include two different blends for the x and y measurements. 
We have also assumed that the blend's proper motion is 
negligible for simplicity, but again it should not be diffi- 
cult to generalize to moving blends. Now, the photometry 
alone gives us t Q , ig, u m i n , F s , and The remaining 6 
parameters are to be derived from the astrometry. First 
note that the unblended astrometric motion along the rel- 
ative proper motion vector is antisymmetric in time with 
respect to ioj while the unblended motion perpendicular 
to it is symmetric (see Figure 6). The blend position and 
brightness are constant, and therefore symmetric. So if 
we multiply a; b s by Ftot , and antisymmetrize in time with 
respect to t , then the blend parameters drop out, and 
the resulting motion is along the unblended proper mo- 
tion vector. So we immediately get (f>. We can estimate 
6*e from £e, u mm and the velocity x at peak. Again, note 
that the blend parameters make a constant contribution 
to the product Ftota?obs, so by taking the difference of this 
quantity evaluated at two different times, the blend pa- 
rameters drop out. We can thus estimate the unlensed 
source position x s by using the estimates for the other 
parameters, and the change in F tot x h s from peak to the 
end of the event. Now we have estimates for all parame- 
ters except the blend position Xb, which we can guess by 
comparing the measured -F to t:Eobs with the contributions 
from the other terms. Thus, we can simply and accurately 
guess values for all of the microlensing parameters, and 
use these as initial estimates for fitting routines. This is 
also true if there are different blend fractions for the two 
baseline orientations. 

To quantify the above statements, we have simulated 
SIM observations of microlensing events. Each of our sim- 
ulated events has f?E = 0.32 mas and £e = 20 days. The 
other parameters, such as u m i n , 4> and x s were picked ran- 
domly. We measure 25 points during the event, evenly 
spaced in time, starting when the apparent magnification 
exceeds a threshold A > A t h = 1.5, and lasting 5£e from 
the trigger. This is probably not an optimal observing 
strategy, but it is not our intention here to optimize the 
strategy. Each astrometric measurement is made following 
the procedure described in previous sections - the fringe 
pattern is simulated and fitted to recover the source po- 
sition and brightness. We assume that the x and y mea- 
surements occur at the same time, although this assump- 
tion has no effect. If SIM measures the orthogonal base- 
lines at different times, say to minimize its overall slewing, 
no information should be lost. Even if the baselines are 
not perfectly orthogonal, measurements of the microlens- 
ing parameters should not be seriously degraded. 

To measure how badly crowding affects microlensing 
measurements, we first see how well the parameters may be 
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recovered for a single, unblended source. To save on com- 
putation time, we used white-light fringe measurements, 
which as noted above should not be much worse than 
spectrally dispersed measurements for unblended sources. 
Now, to measure the lens mass and distance, the param- 
eters of interest are 0e, £e, and <p. In the first panel of 
Figure 7 we plot fit residuals for the above observing strat- 
egy. Note that these fits utilize both SIM astrometric and 
photometric measurements of the microlensing events. We 
find that 6>e can be measured to 2.84%, and ts to 0.5%, 
for input values of 6e = 0.32 mas and £e = 20 days. 

Now, we add blend stars to see how badly they affect 
microlensing mass and distance measurements. Here, we 
clearly need to use the full multifrcquency fringe fit to re- 
solve bright blend stars. From the individual fringe fits, 
which gave errors ~ 10 ^as (twice as large as the unblended 
errors) we might think that the errors here would be twice 
as large. We plot the fit residuals in the second panel of 
Figure 7. We obtain errors a$ = 0.008 mas, a^ = 0.0235, 
and at = 0.142 days. Normalized to 9e = 0.32 mas and 
tE = 20 days, these are fractional errors of 2.5%, 2.4%, 
and 0.71% respectively, as good as or better than the un- 
blended white light errors! Again, the reason that these er- 
rors are not twice as big as the unblended errors is that the 
relatively large single measurement residuals mainly arose 
from unresolved sources falling inside the target's central 
fringe, which as we discussed may easily be removed. 

Another important question to ask is how badly we 
would do without the spectral decomposition of the fringe 
pattern. To address this question, we repeated the above 
simulations, but this time used only white-light fringes. 
The resulting errors are plotted in Figure 7. We find that 
the error distribution is bimodal, and neither blob is cen- 
tered on the correct answer! The bulk of the fits underes- 
timate #e by roughly 10%, with a scatter of ~ ±5%, al- 
though a subset of events give somewhat overestimated O-g, 
about 5% too large, with a scatter of ~ ±3%. This "more 
correct" subset, however, overestimates £e by roughly 8%. 
So if we were to fit a white light measurement of a mi- 
crolensing event in the LMC, we could be reasonably cer- 
tain that the derived masses and distances would be incor- 
rect by about 20%. While 20% distances would certainly 
be good enough to resolve the question of the location of 
the LMC lenses, it is clearly preferable to utilize the chan- 
neled spectrum when possible. 
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Fig. 8. — Fit residuals for microlensing event with an equally 
bright blend in FOV, fitting to the white-light fringe pattern. We 
only plot for distances large compared to the fringe size; closer sep- 
arations give much larger errors, erg is in mas. 



We should be careful, however. Most field blends drawn 
from the LMC luminosity function are less than a few per- 
cent of the brightness of a V — 19 star. A comparably 
bright star landing in the field of view could corrupt mea- 
surements based upon white-light fringes (Rajagopal et al. 
2001). To illustrate this, we simulate white light SIM mea- 
surements of microlensing events with an equally bright 
source in the field of view. Now, it should be clear that 
sources within a few fringes of the target will hopelessly 
corrupt our measurements. However, sources falling out- 
side the coherence envelope have significantly diminished 
fringe contrast, and one might hope that their effect would 
be minimal. These hopes are unfounded, as figure 8 illus- 
trates. This figure plots the errors in #e as a function of 
source separation (both in mas). Even sources lying more 
than 100 mas away can induce 20% errors in the fit. A 
caveat to this statement is that the assumption of rectan- 
gular bandpass overestimates the effect of distant sources. 
The discontinuous jumps in F(k) at the edges of the band- 
pass give sidelobes in the white-light fringe envelope, which 
would be suppressed with a more smooth decline. Never- 
theless, distant bright sources are clearly cause for concern 
when centroiding white light fringes. Now, space telescope 
or ground based AO snapshots could detect such sources. 
However, the positions available from such snapshots are 
rough by SIM standards, ^ 20 — 40 mas, which is larger 
than a fringe size and larger than the length scale of the 
variations exhibited in Figure 8. Thus, such high reso- 
lution snapshots would not be helpful for fitting multiple 
sources. 

How much better do we do if we have dispersed fringes? 
As previous sections showed, bright sources landing out- 
side the central fringe are easily detected and removed. 
Therefore, such sources will not corrupt SIM microlensing 
events beyond adding photon noise, which can be allevi- 
ated by integrating longer. We might worry about closer 
sources, however. Figure 9 shows results from the same 
simulation as above, but for blends within 50 mas, using 
the channeled spectrum. The errors are all small (~ 3%) 
except for marginally resolved blends at about 3 mas sep- 
aration which have 10-20% errors. Bright sources at this 
distance fall within the central fringe and cannot be clearly 
resolved, but are not well enough inside the central fringe 
to be treated as completely unresolved. In contrast, faint 
sources at this distance are effectively unresolved. 

This shows the extent to which spectrally dispersed 
fringes are useful. If the spectral decomposition of the 
fringe pattern is unavailable or costly to transmit, and we 
must rely upon white- light fringes, then high resolution 
space telescope or AO images will be required before fol- 
lowing up a microlensing event with SIM. Detection of a 
bright field source in such snapshots would not help us 
remove the blend's effect, but instead would tell us that 
such events are hopelessly corrupted and should not be 
followed by SIM. Since bright events in the LMC are rare, 
and snapshots would tell us nothing about possible blends 
< 40 mas from the target which would also corrupt SIM 
measurements, we hope that spectrally dispersed fringes 
will be available in the final SIM design. High resolution 
snapshots may be useful with dispersed fringes as well, 
since they could be used to avoid orientations that project 
bright field stars inside the target's central fringe. 
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Fig. 9. — Fit residuals for microlensing event with equally bright 
blend star and 15 other blends drawn from LMC luminosity func- 
tion. The 15 faint LMC blends are placed randomly in field of view, 
and errors are plotted as a function of the bright blend position. 
Compare this curve with Figure 4. 



5. SUMMARY 

We have shown that the spectral decomposition of the 
fringe pattern, which is often used for group delay estima- 
tion, also allows highly precise astrometry to be performed 
in crowded fields. Even when white-light fringes are hope- 
lessly corrupted, dispersed fringes allow microarcsecond 
positions to be measured. We find that any significant 
source outside of the central fringe (> 6 mas) of the target 
star may be detected and removed by fitting the two di- 



mensional (delay vs. wavelength) fringe pattern, thereby 
restoring SIM astrometry to the precision (few ^as) ex- 
pected for isolated sources. We have further shown that 
in mildly crowded stellar fields, precise astrometric imag- 
ing is possible using fringes in the channeled spectrum, 
with positional precision comparable to that expected for 
single star astrometry. Although the details will depend 
upon SIM's precise characteristics, we expect that imaging 
should be possible for up to 10 comparably bright sources 
in the field of view. 

We then illustrated our method with the specific ex- 
ample of astrometric microlensing. We showed that, al- 
though unresolved sources falling within the central fringe 
could corrupt individual astrometric measurements, accu- 
rate masses may still be measured by fitting the entire 
microlensing event. In this manner, confusion or blending 
errors may be brought below photon noise errors. For the 
worst case, i.e. LMC events which are the most crowded 
and have relatively faint sources, 5% mass measurements 
are still easily attainable. We also considered the effect of 
exceptionally bright field stars or binary companions, and 
showed that for all separations except about 3 mas, such 
bright sources also do not adversely affect SIM astrometry. 
Chance separations of 3 mas, however, can lead to 20% er- 
rors in mass and distance measurements from microlensing 
events. 

In conclusion, we find that SIM astrometry will not be 
corrupted by blending, and extremely precise imaging of 
mildly crowded fields may be performed using only two 
orthogonal baseline orientations, allowing microarcsecond 
positional measurements. 
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